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This thesis developed the ballistics processor of a 
multiple processor airborne tactical system. The multiple 
processor system consisted of three INTEL-8080 microcomputers: 
the executive processor, the navigational processor and the 
ballistics processor. The ballistics processor utilized a 
general second order Runge-Kutta method of integration of 
the equations of motion of unguided air- to-surface weapons. 

The ballistics processor computed sufficiently accurate and 
timely solutions to enable the executive processor to extra- 
polate an accurate release point for the weapon. The 
algorithm permitted complete flexibility in release conditions 
and allowed a complete arsenal of air- to-surface weapons 
currently carried on the A7-E aircraft. The cost of the 
ballistics processor using current "off the shelf" components 
is $1635 and the entire tactical system was estimated to cost 
$3891. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 10 

II. BACKGROUND 15 

A. THE NAVIGATIONAL COMPUTER - 16 

B. THE BALLISTICS COMPUTER 18 

C. THE EXECUTIVE COMPUTER 18 

D. THE INTERCOMPUTER COMMUNICATIONS NET - 19 

III. BALLISTICS PROBLEM - 22 

A. BASIC ASSUMPTIONS 22 

B. DERIVATION OF THE EQUATIONS 23 

C. RUNGE-KUTTA APPLICATION 24 

D. CLASSES OF WEAPONS 2 5 

1. Streamlined Bombs and Bullets 25 

2. Drogued Bombs 26 

3. Cluster Bombs 26 

4. Unguided Rockets 26 

E. COMPRESSION OF DRAG DATA 27 

F. CURVE FITTING THE DRAG COEFFICIENTS - 29 

IV. FORTRAN SIMULATION - 38 

A. BOMB DROP SIMULATION 38 

B. OVERALL SIMULATION 4i 

V. PL/M IMPLEMENTATION - 45 

VI. CONCLUSIONS - - 47 



APPENDIX A PROGRAM FLOW CHARTS AND PROCESS GRAPHS -- 
APPENDIX B PROGRAM DEFINITION OF VARIABLES 



69 




_ 



FORTRAN AND SIMULATED PL/M OUTPUT 7 2 

FORTRAN PROGRAM 74 

PL/M PROGRAM 84 

LIST OF REFERENCES - 102 

INITIAL DISTRIBUTION LIST 103 



6 



LIST OF TABLES 



Table 

I. REFERENCE DRAG CURVE COEFFICIENTS AND CUTS 33 

II. WEAPON CONSTANTS IN DECODE 34 

III. SAMPLE ERRORS IN IMPACT RANGE CALCULATIONS 
FOR THE GENERAL-PURPOSE BOMB MK82/SNAKEYE/ 

UNRETARDED - 39 

IV. SAMPLE ERRORS IN IMPACT RANGE CALCULATIONS 
FOR THE SPECIAL-PURPOSE BOMB MK40/MOD 0/ 

CONICAL TAIL 39 

V. SAMPLE ERRORS IN IMPACT RANGE CALCULATIONS 
FOR THE PRACTISE- PURPOSE BOMB MK76/WITH 
LUG - 40 

VI. FORTRAN TEST CASES (TIME INTERVAL, MK83) 42 

VII. PL/M TEST VERIFICATION (TIME INTERVAL, 

MK83) - 46 



7 



LIST OF FIGURES 



Figure 

1. Attack Aircraft's Tactical System's 

Interconnection 11 

2. The Intercomputer Communication Net 20 

3. Curve Fitting Cd 30 

4. Thrust Versus Time for a Typical Rocket 36 



8 



ACKNOWLEDGEMENT 



The author wishes to express his thanks to Associate 
Professor Kodres for his overall assistance and guidance 
throughout the project, and his keen insight into micro- 
computer programming. Special thanks go cut to my family, 
Jeanne and Tim, who have constantly served as an inspiration 
and stabilizing influence throughout the preparation of this 
thesis . 



9 



I. INTRODUCTION 



Military airborne tactical systems are used to aid the 
pilot and bombardier/navigator in their specific functions. 

For an attack aircraft, the tactical system is used in the 
following manner: 

1. to navigate to a given target 

2. to issue steering commands enroute and/or to guide 
the autopilot 

3. to display to the crew the current target location 
in relation to the aircraft 

4. to solve ballistics equations which enable the crew 
to release a weapon at the appropriate instant. 

The presently operational systems employ a small 
general purpose computer, such as the IBM 4 PI system. In 
this thesis, the ballistics processor of a multiple processor 
organization of a typical tactical system for an attack air- 
craft is presented. The system's architecture made use of 
three INTEL- 8080 microcomputers, each of which was dedicated 
full time to its assigned function. The purpose of the 
multiple processor organization is to exploit the presently 
emerging inexpensive microprocessor technology to solve 
military tactical problems. The system's organization is 
shown in Figure 1. Because microcomputers presently have 
no hardware multiply-and-divide functions, the programmed 
software multiply-and-divide functions caused the microcomputers 
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FIGURE 1, ATTACK AIRCRAFT'S 
TACTICAL SYSTEM'S I HTERCO'INECT I Oil 
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to run on a typical mix of operations approximately six times 
slower than the presently operational computer systems. By 
reorganization of the computations and the use of table 
look-up techniques, the described system was capable of meet- 
ing the real-time requirements of an airborne weapon delivery 
system. 

The ballistics computer uses the best estimate of the 
position, dive angle, altitude, and velocity of the aircraft 
in its computation. It is also provided with the weapon 
identification (for the particular drag characteristics) 
and the wind velocities. From this information and the 
target's altitude (above Mean Sea Level), the computer 
determines the down range travel and time of fell of the 
unguided weapon (i.e., the impact coordinates in the target 
plane) . 

The accuracy with which a ballistic weapon can be 
delivered against a target depends greatly upon the accuracy 
of the sensor- suppl ied release parameters as well as the 
accuracy with which the ballistics equations are solved to 
predict the weapon’s down range travel. 

Most aircraft unguided air-to-ground weapons can be 
described as ballistic projectiles. The only forces acting 
on them after release from the aircraft are gravity and 
aerodynamic drag. Bullets, streamlined bombs, drogued 
(retarded) bombs, cluster munitions and unguided rockets 
(after motor burnout) are all ballistic projectiles. Guided 
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weapons and weapons developing lift are not ballistic pro- 
jectiles nor are they considered in this development. 

This thesis first describes the over-all tactical system 
in order to give the reader a perspective of how the ballistics 
processor functions within the entire weapon system. Included 
in this background material is a discussion of the naviga- 
tional computer, the ballistics computer, the executive com- 
puter and the intercomputer communications network. 

The basic calculations involved in solving the ballistics 
problem is described next. The assumptions used in the pro- 
blem formulation are stated; the derivation of the particular 
differential equations is presented; and the application of 
the Runge-Kutta integration method is shown to be suitable 
for current operational airborne weapon systems. The four 
general classes of weapons and how the algorithm handles each 
integration is defined. The unique method of compressing the 
drag data for all twenty- eight weapons is shown to be a simple 
curve fitting of the drag coefficients with respect to seme 
standard. 

The simulation was initially programmed using FORTRAN 
for ease of testing and debugging. The important phases in 
programming were in the DERIV and RUNGE subroutines where a 
close study of the optimum calculation methods was conducted. 
The program was tested for accuracy of results by comparison 
with the published range tables (NAVA I R 01-lC-lT-l). 

Upon completion of the FORTRAN implementation, the re- 
sultant program was translated into PL/M, a high level 
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language compatible with the selected microcomputer 
80). The selected test cases were verified correct 
PL/M version of the program. 



(INTELLEC 
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II. BACKGROUND 



In order to give the reader a complete understanding of 
how the ballistics processor functions, this section is 
devoted to describing the overall tactical system in which 
the ballistics processor is an important part. The first 
section describes the function of the navigational computer, 
the sensors which provide data to the computer and the 
algorithm used to compute present position. 

The ballistics computer will be completely defined in 
the remaining sections of this thesis and will only be 
briefly reviewed here. Using the best estimate of position 
and velocity at a given instant, together with information 
of the weapon’s drag characteristics, the ballistics 
computer solves a system of four differential equations and 
calculates the impact coordinates of the weapon in the target 
plane. 

The executive computer controls the display and generates 
commands to the autopilot and weapon release mechanism. This 
computer is interrupted asynchronously as fresh data are 
developed in the navigational and the ballistics computers. 
The executive computer extrapolates from the data the appro- 
priate release instant of the weapon. 

The intercomputer communication net is described next. 

The star-like communication net. provided a simple and rapid 
means of information transfer. 
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A. THE NAVIGATIONAL COMPUTER 



This subsystem is present in all tactical systems. In 
operational systems, a subprogram (module) is periodically 
executed to update the present position by the change in 
position in the last time increment. 

In this computer, the sensor instruments provide digital 
data to the computer at specified time intervals, controlled 
by a clock. There are four sets of sensor instruments which 
provide the data: 

1. The Inertial Navigation System (INS) 

2. The Doppler Radar (DR) 

3. Air Data (AD) 

4. Radar Set (RS) . 

The INS system provides attitude (roll and elevation) and 
velocity increments in the x, y, z directions. The DR is a 
velocity sensor that uses the doppler principle for the 
continuous measurement of aircraft ground track speed and 
drift angle. The AD computes outputs of modified corrected 
static pressure, pressure altitude and mach number from 
static and pitot tube pressure inputs. The RS provides 
target azimuth, elevation, and range signals to the naviga- 
tional computer. A magnetic compass measures the aircraft 
heading in relation to magnetic north. In order to navigate 
accurately in this mode the wind speed and direction are 
determined by using the "sensed" velocities of the aircraft. 
Altimeter measurements are also sampled in this system to 
update the present position altitude. 



16 



There are four basic modes in which this system operates. 
The crew has a choice of using both the INS (velocities) and 
DR outputs to minimize the Schuler oscillations. This is the 
most reliable and accurate mode. If the crew suspects that 
either the INS or the DR system is malfunctioning, they may 
use only the data from one of the sources. The least 
reliable and inaccurate mode is the last resort, when both 
the DR and the INS are malfunctioning. The AD mode uses the 
externally estimated wind vector and true airspeed of the air- 
craft to determine the position in latitude and longitude. 

In all operating modes, the same basic technique of smooth- 
ing input data is used. The least squares curve fitting 
technique is used to fit a quadratic polynomial to the most 
recent n (n ^ 16) data values. Use of the Legendre 
orthogonal polynomial functions is made in order to make the 
numerical calculations as accurate and rapid as possible. 
Smoothing the short time increment data values will generate 
less total error in the results. 

Using the smoothed data values, the appropriate calcu- 
lations are made to estimate the distance increments in the 
past time interval. As soon as the new value is generated, 
the display of the navigational computer is refreshed, the 
executive computer is interrupted and this most recent 
estimate is transferred to the executive computer. If the 
navigational computer is functional, the crew has information 
on the present position, independent of the executive 
computer. 
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The least squares curve fitting technique enables the 
system to diagnose sensor malfunction quite easily. If the 
sensor generates fluctuations in the input data above a 
certain threshold, then the sensor can be diagnosed as mal- 
functioning. Similarly, if two sensors disagree above an 
allowable threshold, then an automatic data analysis will 
suggest which sensor is failing. If the crew concurs in 
this diagnosis, the operating mode of the navigational 
system can be changed to exclude the faulty sensor. 

B. THE BALLISTICS COMPUTER 

The ballistics computer is provided with the present 
estimates of the position, dive angle, altitude, and 
velocity of the aircraft. It is also provided with the type 
of weapon and wind conditions. From this information and 
the target's altitude, the computer determines the down- 
range travel of the weapon, if released at this instant. 

The assumption of no guidance after release from the air- 
craft is still applicable. However, the weapon may have a 
rocket assist, or a retarded mode of fall. 

C. THE EXECUTIVE COMPUTER 

The executive computer controls the main display and 
generates steering commands to the autopilot and firing 
pulses to the weapons release mechanisms. One of its 
critical functions is to extrapolate from the data generated 
by the ballistics computer the weapon release point. 
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A real time clock is used to record the time at which 



the asynchronous data were generated by the ballistics 
computer. The position of the impact point in the target 
plane can be described by a quadratic function of time. 
Therefore, the real time release point (adjusted for time 
delays) can be extrapolated very accurately from these data 
without additional extensive calculations. 

D. THE INTERCOMPUTER COMMUNICATIONS NET 

In order to keep the design of the intercomputer 
communications net simple and inexpensive, the starlike 
network was chosen. The executive computer is the center of 
the star and one input port is multiplexed to the two output 
ports of the navigational and ballistics computer as shown 
in Figure 2. 

The intercommunication takes place in the following 
manner: An initializing phase of the program interrogates 

the keyboard to obtain data about the initial position of 
the aircraft. The latitude and longitude are entered in 
degrees and minutes to the nearest 1/10 minute; the 
altitude is initialized at the height above mean sea level. 
These data are made available to the navigational computer. 
Also the type of weapon and the position of the target (if 
known in advance) is made available to the ballistics 
computer. Similarly, the mode of navigation is determined 
(if possible prior to flight) and relayed to the navigational 
computer . 
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FIGURE 2, THE INTERCOMPUTER 
COMMUNICATION NET 
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The execution phase of the computer program makes use of 
the real time clock in the executive computer which is 
asynchronously interruping the navigational computer in 
order that sensor information be transferred via input port 2 
to the navigational computer. When the navigational computer 
has completed its calculations, it interrupts the executive 
computer and transfers the calculated current position 
sequentially via input port 1 to the main computer. The 
amount of information sent is 9 bytes, and the transfer is 
accomplished in about 250 micro- seconds . 

Similarly, when the ballistics computer has completed 
its calculations, it interrupts the executive computer and 
transfers its information sequentially via input port 1 to 
the executive computer. Only 6 bytes of information are 
sent in approximately 170 micro- seconds . Between interrup- 
tions, the executive computer updates the display using 
output port 2. Calculations are continuously made to extra- 
polate the release instant of the weapon, and the release 
mechanism is activated when appropriate via output port 3. 



III. BALLISTICS PROBLEM 



The equations of motion governing the trajectory of a 
ballistic projectile are a simple- appearing set of second- 
order differential equations whenever the drag is neglected 
However, when a reasonably accurate model of the aerodynami 
drag caused by the projectile's motion through the air mass 
is included in these equations they are rendered nonlinear, 
and no general solution in closed form has been found that 
is satisfactory for solution by an airborne fire-control 
computer. 

A. BASIC ASSUMPTIONS 

The mathematical model was chosen to represent reality 
as accurately as possible. The armed forces publish range 
tables for various weapons and therefore the same model was 
chosen as the one used to construct these tables (NAVAIR 01 
1C- IT- 1) . 

The equations of motion were developed assuming the 
projectile is a point mass acted on only by the force of 
gravity and the retardation forces due to air resistance. 
The trajectory was restricted to a plane by ignoring the 
crosstrail effects of the wind. These crosstrail effects 
and the component of wind which affects drag are considered 
in the executive computer. 
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The assumptions are summarized below: 

1. The Earth is flat and nonrotating. 

2. The gravitational attraction is constant. 

3. The projectile is a point mass. 

4. The projectile is not powered (except rockets) and 
has a constant mass. 

t 

B. DERIVATION OF THE EQUATIONS 

The weapon’s fall is assumed to be governed by a drag 
force which is proportional to the square of the velocity. 
Newton's equation of motion becomes 

m dVt/dt = - g m - || Vt|| 2 Cw Vt/ ||vt|| 

where the only two forces acting on the weapon are gravity 
and aerodynamic drag. 

Vt - total velocity 
m - mass of the projectile 
g - force due to the Earth's gravitation 
Cw - drag coefficient which depends on the 

characteristics of the weapon and may vary 
with time (in case the weapon has two drag 
curves), and it varies slightly with airspeed. 
In order to solve this differential equation, initial 
conditions, that is the position and velocity, must be 
specified. The most convenient coordinate system is height 
"y" above sea level and "x" downrange travel, which is zero 
at the time of release. 
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In the x and y component form, Newton's equation becomes 



dVx/dt 

dVy/dt 

dx/dt 

Dy/dt 

x(0) 

y (o) 

Vs (0) 
Vy(0) 



- || Vt || /m * Vx * Cw 

- g - || Vt || /m * Vy * Cw 
Vx 



Vy 

0 



downrange travel 
present altitude 
present horizontal velocity 
present vertical velocity 



C. RUNGE-KUTTA APPLICATION 



The choice of the "best" Runge-Kutta (R-K) formula for 
a particular problem was not entirely obvious. The lower 
order formulae are simple, but require more integration 
steps than the more accurate higher order formulae. On the 
other hand, the complexity of the higher order formulae 
may require so much computer time that this becomes the 
overriding consideration. Generally, it was found to be 
safer to use the lower order formulae when discontinuities 
in the coefficients of the differential equations can be 
expected. Such discontinuities frequently occur in bomb 
trajectories (e.g., drogue deployment or thrust termination). 
The second-order R-K formula also afforded the greatest 
flexibility in the distribution and size of the integration 
steps. These considerations, along with their simplicity, 
led to the choice of the second-order R-K method as the most 
appropriate for the ballistics algorithm. 
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Let U1 



Vx 



FI = -1/m 



* Vx * Cw 



U2 = y 
U3 = Vy 




F2 = Vy 

F4 = -1/m * ||vt|J * Vy - Cw 



g 



The differential equations may be expressed as 

dUi/dt = Fi(Ul,U2,U3,U4) i = 1,2, 3, 4 



Therefore, the Runge-Kutta second order method used in the 
program becomes 

mli = At * Fi (U1 ,U2 ,U3 ,U4) i = 1,2, 3, 4 

m2i = At * Fi(Ul + cAtmll,U2 + cAtml2,lJ3 + cAtml3, 

U4 + cAtml4) c = 0.7, i = 1,2, 3, 4 

U’i = Ui + 1/ 2c * ((2c - 1) * mli + m2i) i = 1,2, 3, 4 



The procedure computes new values for the variables 
until the remaining time of fall, DTV, (as computed for a 
vacuum release) , is less than the value of the next integra- 
tion time increment. The procedure then uses the value of 
DTV for the final integration step size, and in the calcula- 
tion of the down range travel of the weapon. 

D. CLASSES OF WEAPONS 

1 . Streamlined Bombs and Bullets 

The process described above represents streamlined 
bombs quite well. After being fired, a bullet is really just 
a small, streamlined bomb and was handled as such. The muzzle 
velocity (3300 feet per second) of a particular gun-bullet 
combination was added to the aircraft velocity. 
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2 • Drogued Bombs 

These bombs are more complex in their trajectory 
modeling. Their common characteristic is that they are all 
released in a relatively low-drag configuration. At some 
time after release, they deploy vanes or a parachute which 
greatly increases the total drag. To compensate for this 
change of drag characteristics, the initial integration step 
size was chosen to coincide with the deployment of the high 
drag devices. Thereafter, there was no difficulty in calcu- 
lating the trajectory of the drogued weapons. 

3 . Cluster Bombs 

These weapons are released as large, low-drag 
containers. At a predetermined point in their trajectory, 
which is time dependent, they use some method to dispense 
smaller weapons. These smaller weapons are of higher drag 
than the container. The ballistics computer will predict 
the impact point of the center of the pattern rather than 
individual positions of the bomblets. The method of applying 
the algorithm is to calculate the trajectory of the container 
to the point where it dispenses the small weapons; switch 
to a different drag curve which describes the motion of the 
pattern center; and integrate down to the target altitude. 

As in the retarded bombs, the initial integration step was 
chosen to coincide with the act of dispensing the small 
weapons . 

4 . Unguided Rockets 

Rockets present several unique complexities. They 
change weight during flight; they have thrust as well as 
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drag; and they slew around after firing, because their 
launches may not line up with the aircraft direction of 
flight. The algorithm treats unguided rockets like a 
cluster bomb with the first stage having an average thrust 
as well as an average thrust time. This provides sufficient 
accuracy and saves computer storage that would be taken 
storing the thrust profile of each rocket. 

E. COMPRESSION OF DRAG DATA 

The behavior of a particular weapon's drag coefficient, 
Cd, as a function of Mach number must be known to calculate 
impact ranges for that weapon. Normally, this information 
is given in tabular form with some weapons having more than 
100 points in the table. Most aircraft have the capability 
to carry many different types of weapons. It is desirable 
to store all the Cd information of an aircraft's weapon 
repertoire in the airborne computer to avoid loading the Cd 
data for each different combination of weapons. This 
requires storing drag coefficient data in a more efficient 
form than a table. 

The following factors influence the choice of a scheme 
to approximate the drag coefficients for a given set of 
weapons : 

1. Computer storage 

2. Computer computation time 

3. Delivery envelope of each weapon 

4. Allowable downrange and time-of-fall errors. 
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It would be ideal to use one general, fairly powerful 
expression to fit all the drag tables for the weapon 
repertoire of a given aircraft. This is quite difficult to 
accomplish because weapons like guns and rockets require 
Cd values for Mach numbers much higher than free-fall bombs. 
As a consequence of this generally acceptable tabling pro- 
cedure, extra storage is spent on weapons that do not require 
the extra capability. 

A significant saving can be realized when it is 
recognized that several free-fall weapons have drag 
coefficient curves that differ only by a multiplicative 
factor. The Mk SO low-drag series are one such group of 
weapons with the Mk 84 drag curve conventially taken as the 
reference. Similarly, the Mk 106, CBU bomblet, and Sadeye 
bomblet are another group with the Garve drag curve as the 
reference. The specific region of interest of many more 
drag curves can be approximated by multiplying these two 
reference drag curves by a factor and then translating them 
along the Cd and M axes until they match the original drag 
curves . 

To further aid in the explanation of approximating 
weapon drag curves by this method, let Cd(M) and M be the 
drag coefficient and Mach number of a weapon drag curve. 

Let £d (M) be the drag coefficient of a reference drag curve 
and let A, B and DM be the multiplicative factor, drag 
coefficient translation, and Mach number translation , respec- 
tively, that when applied to the reference drag curve will 
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approximate the particular weapon drag curve of interest. 
That is: 



Cd (M) ± A * Cd(M + DM) + B 

A hypothetical weapon drag curve, Cd, and reference drag 
curve, Cd, are shown in Figure 3. For simplicity of discus- 
sion, the Mach number intervals (M2 -Ml) and (M4 - M3) are 
equal. Suppose the Cd curve is multiplied by an A so that 
the Cd and the A * £d curves are nearly identical. For the 
Cd and A * Cd curves to be a congruent, the A * C d curve must 
be translated along the Mach axis an amount equal to DM, and 
along the Drag Coefficient axis an amount equal to B. 

For many weapons, the above stated equation will not give 
an accurate fit over a wide interval of Mach numbers; instead, 
A, B and DM should be chosen so that the Cd fit is best for 
the Mach number interval that corresponds to the delivery en- 
velope for a particular aircraft. 

Considering the weapons currently implemented on a typical 
attack aircraft, (A7-E), both the Mk 84 and the Garve reference 
drag curves are necessary to handle the free-fall weapons. 
Instead of extending these curves to handle guns and rockets 
that have a much larger Mach number range, a rocket-gun 
reference drag curve is used. 

F. CURVE FITTING THE DRAG COEFFICIENTS 

In the previous section, the reference curves needed and 
how they were used to approximate weapon drag curves for use 
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DRAG COEFFICIENT 
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in the algorithm was presented. In this section the reference 
drag curves versus Mach number is represented. The form of 
the function used will minimize computer storage, logic and 
computation time for the allowable error in downrange and 
time of fall for the different weapons. It was decided that 
each of the three reference curves be divided into three 
regions with Cd expressed as a second-order polynomial in 
Mach number for each region. 

Cd CM) = aO + al * M + a2 * M 2 

The coefficients, aO, al, and a2 depend on the region the 
Mach number M is in. The drag coefficient for a particular 
weapon is expressed as: 

Cd (M) - A * (aO + al * (M + DM) + a2 * (M + DM) 2 ) + B 

As a further step to save computer storage and reduce the 

computation time for the algorithm, another compression of 

drag data was undertaken. The computation C * (PI/8 * Cd) , 

CKDG, is made in the Runge-Kutta portion of the algorithm 

2 

where C is the bomb factor, d /W, where d is the weapon 
diameter in feet and W is the mass in slugs. One constant 
per weapon in computer storage and two multiplications per 
integration step is saved by expressing CKDG as 

CKDG = CF * CKDG + DKG 

where CF = A * (d/d) 2 * W/W 
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CKDG = dVW * PI/8 * Cd 
DKG = d 2 /W * PI/8 * B 

In the above equations, d and W are the assumed diameter 
and mass for the reference curve Cd. The final form of CKDG 
in the algorithm is 

CKDG (M) = CF * (bO + bl * (M + DM) + b2 * (M + DM) 2 ) + DKG 

Table 1 gives a list of bO, bl, and b2 for each region of the 
three reference curves and the values of CF, DM, and DKG, are 
tabulated in Table 2. 

Thrust appears in the algorithm in the total drag 
function 



IIH = TH/V - RHO * CKDG * V 
where TH = thrust/mass 
RHO = air density 
V = weapon velocity 

and CKDG was previously defined in the above equations. 
Figure 4 shows thrust versus time for a typical rocket 
approximated by a constant, thrust = total impulse/Tl, 
where T1 is decreased from actual motor burn time in order 
to make thrust equal to the average thrust of region 2. The 
mass of the rocket, W, obviously decreases during rocket 
motor burn time. However, it is adequate for W to be 
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TABLE 1. Reference Drag Curve Coefficients and Cuts. 



Coefficient 


Reference Curve 


Mk 34 


G a r ve 


Ro eke t -Gun 


Region 1 : 










bO 




1 . 5729E-03 


3 .535039 


0 . 104 115 


bl 




0 . 0 


-3.347782 


-0 .230347 


b 2 




0 . 0 


2 . 872624 


0 . 167644 


Region 2 : 










bO 




4 . 67840E-02 


11.26165 


-0. 194037 


bl 




-0.109711 


-27.41625 


0.401478 


b 2 




6 . 65480E-02 


21 . 73083 


-0 . 1646 12 


Region 3: 










bO 




-0 . 11638015 


-23 . 79 154 


7 . 33246E-02 


b 1 




0 . 21764389 


44 . 26077 


-2 . 03275E-02 


b 2 




-9 . 76706E-02 


-14.49960 


2 . 446 82E-03 


Cut 








First cut, 


CTl 


0.834 


0.622 


1 .032 


Second cut , 


CT2 


0.977 


0 .885 


1 . 300 



TABLE 2. Weapon Constants in DECODE 
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TABLE 2. Weapon Constants in. DECODE (continued) 
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FIGURE L \. Thrust Versus Time Curve for a Typical Rocket 



approximated by a constant, by varying W 
error is minimized. The second stage of 
after weapon release greater than Tl, is 
will be zero. 



until downrange 
the rocket, for time 
free fall and TH 
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IV. FORTRAN SIMULATION 



A. BOMB DROP SIMULATION 

The initial ballistic trajectory algorithm for digital 
airborne fire control computers used in the simulation of 
Traj , Deriv, and the Runge routines, was developed by the 
Boeing Company, Seattle, Washington. The research was 
completed under a Naval Weapons Center contract in September 
1970. 

As previously mentioned, the Second-Order Runge-Kutta 
integration scheme was chosen because of its flexibility 
and simplicity. The results, as evidenced in Tables 3, 4 
and 5, compared favorably to the published data (NAVAIR 
01-lC-lT-l), with mean errors of 0.10% for the down range 
travel. The low-drag bombs shewed relatively little 
sensitivity to the step size, so long as the maximum 
integration step, (DMAX) , was less than or equal to 5 
seconds and the normal Initial Integration Step, (DTI), 
was less than or equal to 3 seconds. However, the dual 
stage weapons (IDNO>17) did show a large sensitivity to 
the step size used during the integration scheme. The 
accuracy desired, requii'ed the DTI value to be consistently 
less than or equal to 2 seconds with the average value of 
1.05 seconds, for the final eleven weapons. 

The ultimate verification of the step size for the 
Runge-Kutta (Deriv also) method was conducted by varying 
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TABLE 3. Sample Errors in Impact Range Calculations for the 
General-Purpose Bomb Mk 82 / Snakeye /Unr e tarded . 



Release 


Conditions 


Impact Ranee . f t 


Impact. Range 
error, 
ft 


Velocity , < 

kno t s 


ilt itude , 
ft 


Angle , 
deg 


Ballistic 

Tables 


A 1 g o r i t hm 


400 


2500 


0 


8039 


8039 


0 




5000 


-10 


9158 


9 159 


+ 1 




12000 


-30 


10536 


10537 


+ 1 




12000 


-45 


7530 


7531 


+ 1 


450 


2500 


0 


8998 


8995 


-3 




5000 


-10 


10025 


10023 


-2 




12000 


-30 


11340 


1134 1 


+ 1 




12000 


-45 


7998 


7999 


+ 1 


500 


2500 


0 


9946 


9940 


-6 




5000 


-10 


10843 


10839 


-4 




12000 


-30 


12067 


1206 7 


0 




12000 


-45 


8406 


8407 


+ 1 



TABLE 4. Sample Errors in Impact Range Calculations for the 
Special-Purpose Bomb Mk 40/Mcd 0/Conical Tail. 





<3ond i t ioi 


13 


Imn act Ranee, ft 


Impact Range 
error , 
ft 


Velocity, 
knots 


\lt i tude , 
ft 


Angle , 
dee 


Ballistic 

Tables 


Al g o r i t hm 


400 


2500 


0 


8305 


8301 


-4 




5000 


-10 


9434 


9432 


-2 




12000 


-30 


10847 


10846 


-1 




1 2000 


-45 


769 1 


7690 


- 1 


450 


2500 


0 


9329 


9324 


-5 




5000 


-10 


10342 


10339 


-3 




12000 


-30 


11670 


11669 


-1 




12000 


-45 


8160 


8159 


-1 


500 


2500 


0 


10348 


10342 


-6 




5000 


-10 


11200 


11197 


-3 




12000 


-30 


12401 


12402 


+ 1 




12000 


-45 


8559 


8559 


0 
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TABLE 5. Sample Errors in Impact Range Calculations for the 
Pract ise-Purpose Bomb Mk 76/ with Lug. 



Release Conditions 


Impact Ranee, ft 


Impact Range 
error 

ft 


Velocity , 
knot s 


Vltitude , 
ft 


Angle , 
deg 


Ballistic 

Tables 


Algo r i thm 


400 


2500 


0 


7633 


76 17 


-16 




5000 


-10 


8719 


8710 


-9 




12000 


-30 


9991 


9992 


+ 1 




12000 


-45 


7232 


7235 


+ 3 


450 


2500 


0 


8463 


8459 


-4 




5000 


-10 


9493 


9487 


-6 




12000 


-30 


10712 


10720 


+8 




12000 


-45 


7669 


7674 


+5 


500 


2500 


0 


9268 


9255 


-13 




5000 


-10 


10195 


10192 


-3 




12000 


-30 


11337 


11340 


+ 3 




12000 


■ -45 


8038 


8039 


+ 1 
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the time interval from one to five seconds. The TRAJ 
integration scheme was simultaneously tested using several 
release conditions to demonstrate the model's flexibility in 
simulating bomb loft, level flight release and dive bombing 
conditions. The emphasis, at this point, was on the number 
of functional calls to RUNGE for the projectile to reach the 
target altitude. The results, by altering At, showed changes 
ranging from 0.00% to 0.44% for the time of fall, and 0.00% 
to 0.46% for the downrange travel, with a substantial decrease 
in the number of calls to RUNGE also evidenced (Table 6) . 

B. OVERALL SIMULATION 

The PRINT routine was a means of simulating the various 
cockpit displays in the tactical attack aircraft. The output- 
consisted of the release altitude, dive angle, and airspeed, 
as well as the resulting downrange travel and the time of 
fall of the weapon. 

The INPUT routine simulated the collecting and passing 
of the sensor data to the various processors. The sensors 
included in the attack aircraft are the Inertial Navigation 
System, Doppler Radar, Air Data and the Radar Set. Specifi- 
cally, the velocity, dive angle, altitude and the identifi- 
cation number of the weapon (usually provided via the armament 
panel) are required for the execution of the program. 

The SETDAT routine completed the interface between the 
executive computer and the ballistics computer. Its purpose 
was that of assigning the constants their specific nonvarying 
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TABLE 6. Fortran Test Cases (time interval t mk 33) 
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TABLE 6. Fortran Test Cases (continued) 
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values and reinitializing the variables to zero at the start 
of each new weapon's integration. The routine also performed 
various calculations including the determination of the 
velocity in the x and y directions, the conversion of raw 
velocity to feet per second and the calculation of the release 
angle in radians. 

The weapon dependent constants required for the algorithm 
were stored in DECODE. Also, DECODE contains the necessary 
logic to load the correct constants for the particular 
weapon the crew has selected. After the first pass, there 
is no need to go through DECODE before entering the algorithm 
for a trajectory calculation. The actual scheme used for 
DECODE in reality will depend on the particular aircraft and 
its fire control requirements and how the algorithm is inter- 
faced in the airborne computer. 
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V. PL/M IMPLEMENTATION 



To implement the model on the INTEL 8080 microcomputer, 
the floating point arithmetic package was used to maintain a 
four hexadecimal digit accuracy during calculations. The 
floating point variable therefore occupies three eight-bit 
words: two words for the variable mantissa and one word for 
the mantissa sign bit and the variable exponent. The man- 
tissa of the variable is viewed as a normalized binary number 
with the left most digit always one. 

In order to manipulate these variables, an arithmetic 
subroutine library was supplied by Professor Kodres. At 
present, there exist routines that add, subtract, multiply, 
divide and take the square roots of these floating point 
variables . 

The actual PL/M version of the ballistics program was 
translated directly from the working FORTRAN version. All 
the associated flow charts, process graphs, program listings, 
and the simulated output can be found in the appendix. 

The PL/M routines have been compiled and tested. The 
results compare very favorably with the FORTRAN results. 

The PL/M verification of the test parameters is included in 
Table 7. 
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TABLE 7. PL/M Test Verification (time interval, MK 83) 



RELEASE 

CONDITIONS 


BALLISTIC 

TABLES 


DT = 3.0 


COMPUTATION TIME 


Trial 1 


Tri al 2 


10° 

350 kts 
1000 ft 


11.71 

6701. 


11.69 -. 11 % 
(4) 

6694. 5-. 09% 


1.2 

(sec) 


1.1 

(sec) 


10° 

350 kts 
3000 ft 


17.29 
9826 . 


17.28 -.05% 
(6) 

9820. 5-. 05% 


1 . 8 


1.7 


10° 

650 kts 
500 ft 


13. 89 
14192. 


13.82 -.49% 
(5) 

14140. - . 36% 


1.3 


1.2 


10° 

650 kts 
3000 ft 


20 . 85 
20859 . 


20.81 - . 19% 
(7) 

20820. - . 18% 


2.2 


2 . 1 


0° 

300 kts 
3500 ft 


14 . 85 
7382. 


14.83 -.13% 
(5) 

7379. -.03% 


1.6 


1.4 


0° 

300 kts 
15000 ft 


30.98 
15164 . 


30.95 -.09% 
(11) 

15156 . - .04% 


3 . 2 


3 . 0 


0° 

450 kts 
5000 ft 


17.82 

13127. 


17.81 -.05% 
(6) 

13122 . - .03% 


1.8 


1.7 


-20° 

400 kts 
5500 ft 


12.81 

7963. 


12.80 -.07% 
(5) 

7963. 0.00% 


1.2 


1.1 


-20° 

550 kts 
7000 ft 


13.50 

11387. 


13.47 -.22% 
(5) 

11388. +.01% 


1 . 2 


1.2 


-40° 

450 kts 
5000 ft 


8.18 

4689. 


8.18 0.00% 
(5) 

4689. 0.00% 


0.9 


1 . 0 


-60° 

400 kts 
7000 ft i 


9.61 

3194. 


9.60 -.10% 
(4) 

3195. +.02% 


_ 

0 . 8 


0. 8 
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VI. CONCLUSIONS 



The objective of this thesis was to demonstrate a micro- 
computer based digital system which will perform the same 
functions as the attack portion of a presently implemented 
airborne tactical system. The main purpose of the project 
was to verify the feasibility of implementing the most com- 
plex arithmetic task of the airborne tactical computer by 
use of an INTEL 8080 microcomputer at an estimated cost of 
$3891 for the entire tactical system. 

Because of the rapidly decreasing hard;vare costs of the 
microcomputer and significant increases in processing speeds 
which the microcomputer industry is presently realizing, 
these new, innovative systems, which are in all ways equivalent 
to the presently available avionics systems, can be developed 
at a hardware cost which is nearly 10 times smaller. The 
software development costs are also likely to be smaller 
because the system can be subdivided into smaller, independ- 
ently operating subsystems, each of which will perform a well 
defined task. 

The lack of hardware multiply- and-divide operations in 
the microcomputer is compensated by careful process analysis 
to reduce the number of multiply-and-divide operations. This 
allows three "slow" (in the computer sense) computers to 
solve the same pioblem that one "fast" computer is presently 
solving, at a great hardware cost reduction. The necessity 
for careful analysis of the problem may even disappear soon, 
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because the microcomputer manufacturers are likely to 
incorporate the hardware multiply-and-divide operations. 

This thesis demonstrates the feasibility of using the micro- 
computer as a systems’ building block for airborne tactical 
systems . 



APPENDIX A PROGRAM FLOW CHARTS AND PROCESS GRAPHS 



Symbols Used In Flow Charts 

CZZZ) 



Start/Stop 




Output 




Card 



Input 




O 



Processing Symbol 
(Arithmetic Operation) 



Conditional Symbol 



Subroutine Call 




Return 



Statement 



o 
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Program Entry Point (label 



Symbols Used in Process Graphs 

input 

1. Bite Operations (8 bits): 

ADD 

MULTIPLY 
SUBTRACT 
DIVIDE 
ASSIGNMENT 

2. Floating Point Operations: 

ADD (Subtract) Routine 
MULTIPLY Routine 
DIVIDE Routine 
SQUARE ROOT Routine 

3. Subroutine Calls 



ADD 



MULT 



DIV 



5QRT 




)ERIV 



First Derivative Routine 



output 
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MAIN PROGRAM 
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MAIN PROGRAM (PART 2) 




2 



THETA=DEG ( I ANG) *RAD 
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MAIN PROGRAM (PART 3) 



2 




-o© 

-o© 

-£>© 

-o© 
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INPUT SUBROUTINE 






54 • 



PRINT SUBROUTINE 
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SETDAT SUBROUTINE 





CFORMl= 0.0 
CFORM2= 0.0 
DM1= 0.0 
DM2= 0.0 
DKG1= 0.0 
DKG2= 0.0 
VMUZ= 0.0 
VE= 0.0 
SL= 0.0 
DMAX= 5.0 
*** ITYPE= -1 
FN= 0.0 
I BOTH = 1 
* *T= 0.0 

IE 



U= VKTS(IV) *1.6878 
DEL= A TAN (VE , U) 

V=SQRT (U*U+VE*VE) 

VXA= (V+VMUZ) *COS (THETA-DEL) 
VYZ= (V+VMUZ) *SIN (THETA-DEL) 




RETURN 



* SET VALUES FOR PL/M ARE 1,2,3 
** T IS TM AND A IS A IN PL/M. 

*** ITYPE HAS A VALUE 3 FOR -1 IN PL/M. 



I 



56 



DECODE SUBROUTINE 
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DECODE SUBROUTINE (PART 2) 



o© 






E> 



SB 



DECODE SUBROUTINE (PART 3) 





I REF= 1 
I TYPE— 1 
IBOTH-' 2 
DKG1= 7.33E-03 
CFORM2= 1.69E-02 
DM2= 0.38 
DKG2= 0.17166 
DS= 0.6617 
SL= -0.000269 
DTI= 2. 



IREF= 1 
ITYPE= 1 
I BOTH = 2 
CFORMl= 2.0754 
CFORM2= 0.2217 
DS= 4.267 
DTI= 1.5 
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DECODE SUBROUTINE (PART 4) 




60 





61 




3 IN PLM 



62 





63 



RUNGE SUBROUTINE 



* A 

* * <p 



IS TM IN 
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RUNGE PROCESS CHART 



£-• 

§ 




o 

r-H 

o 
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runge process chart 



£ 

< 

04 




o 



r- oo 
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DERIV SUBROUTINE 




V= SQRT ( VX*VX+VY *VY ) 

CM- V* (8.9E-04+3,2E-09*Y)+DM 




1 c> \7 <> 

CKDG= DKG+CF* ( CC ( I REG , 1 , MSTG) + ( CC ( I REG , 2 , MSTG) +CC ( I REG , 3 , MSTG) *CM) *CM) 

V ~ 

HH= TH/V-RHO*CKDG*V 
AN2= HH*VX 
AP2= HH *VY-G 
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DERIV PROCESS CHART 
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MULT 



APPENDIX B: PROGRAM DEFINITION OF VARIABLES 



SYMBOL 

PL/M if 



FORTRAN 


varied 


DEFINITION 


A 


A1 


Runge Kutta parameter 


AA 




0.5/Runge Kutta parameter 


ALT 




Weapon release altitude 


AN 1 , AN 2 




Runge Kutta variables (x direction) 


API } AP2 




Runge Kutta variables (y direction) 


CC 




Matrix of drag curve coefficients 


CF 




CKDG stretch factor 


CFORM1 




CKDG stretch factor for the first stage 


CFORM2 




CKDG stretch factor for the second stage 


CKDG 




Bomb coefficient times drag coefficient 


CM 




Mach number 


CT 




Matrix of Mach curve coefficients 


DEG 




Weapon release angle 


DEL 




ArcTan (ejection velocity/weapon velocity) 


DKG 




Shift in CKDG 


DKG1 




Shift in CKDG for the first stage 


DKG2 




Shift in CKDG for the second stage 


DM 




Shift in Mach number 


DMAX 




Largest integration step size 


DM1 




Shift in first stage Mach number 


DM2 




Shift in second stage Mach number 


DS 




Integration step size 
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SYMBOL 

PL/M if 

FORTRAN varied DEFINITION 



DTI 


Initial integration step size 


FN 


Equal to the Thrust for the first stage 


FRACT 


Parameter for Runte Kutta integration 
(value between zero and one) 


G 


Acceleration due to gravity 


HH 


Total drag function 


IANG 


Number of release angles integrated 


I BOTH 


Key determining if two drag curves 
are required for the weapon 


IDNO 


Identification number of each weapon 


IREF 


Reference for Mach curve coefficients 


IREG 


Mach region index 


IV 


Number of velocities integrated 


IY 


Number of altitudes integrated 


MSTG 


Weapon stage index 


RAD 


Conversion factor for radians 


RHO 


Air density 


SET 


Key to initializing variables 


SL 


Slope coefficient for retarded snakeye 
fins deployment time 


T TM 


Elapsed time of flight from weapon release 


TH 


Thrust/Mass of rocket 


THETA 


Adjusted weapon release angle 


U 


Aircraft speed 


V 


Velocity of weapon 


vx 


Velocity component of V in the x 
direction 
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SYMBOL 

PL/M if 

FORTRAN varied 

VXA 

VXO 

VY 

VYA 

VYO 

VZ 

X 

Y 

YO 

YT 



DEFINITION 

VX at weapon release 

VX at start of present integration step 

Velocity component of V in the y 
direction 

VY at weapon release 

VY at start of present integration step 
Muzzle velocity of guns 
Weapon ground range from release 
Weapon altitude above sea level 
Y at start of present integration step 
Target altitude above sea level 
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FORTRAN AND PL/M RESULTS 



I DNO 


VKTS 


ALT 


DEG 


FORI 


[RAN 


PL/M 


TOF 


RANGE 


TOF 


RANGE 


1 


450. 


3000 . 


-10. 


10.26 


7515. 


10 .27 


7518 . 


2 


450 . 


3000 . 


-10. 


10.41 


7393. 


10.41 


7395 . 


3 


450. 


3000. 


-10. 


10.32 


7467 . 


10 . 32 


7470. 


4 


450. 


3000. 


-10. 


10 . 73 


7189. 


10.74 


7191. 


5 


450. 


3000. 


-10. 


10.91 


7085 . 


10 .87 


7078. 


6 


450. 


3000. 


-10. 


11 . 13 


6955 . 


11. 14 


6957. 


7 


450 . 


3000. 


-10. 


10.32 


7466. 


10.33 


7469. 


8 


450. 


3000. 


-10. 


10.59 


7292. 


10 . 60 


7294 . 


9 


450 . 


3000 . 


-10. 


10.29 


7492. 


10 . 29 


7495 . 


10 


450. 


5000 . 


-10. 


10.25 


7521. 


10 . 26 


7524. 


11 


450 . 


3000. 


-10. 


10.25 


7529. 


10.26 


7552. 


12 


450. 


3000. 


-10. 


10.24 


7536. 


10 . 24 


7539 . 


13 


450 . 


3000. 


-10. 


10 .22 


7547. 


10 . 23 


7550. 


14 


450. 


3000 . 


-10. 


10.31 


7478. 


10 , 31 


7480. 


15 


450. 


3000 . 


-10. 


10.40 


7424. 


10.41 


7427. 


16 


450. 


3000. 


-10. 


10.26 


7516. 


10.27 


7 519. 


17 


450. 


3000. 


-10. 


10.50 


7368. 


10 .51 


7571. 


18 


450 . 


3000. 


-10. 


16.55 


3852. 


16.57 


3847. 


19 


450. 


3000. 


-10. 


20.84 


4620. 


20.87 


4619. 


20 


450 . 


3000. 


-10. 


14.92 


5771 . 


14.94 


5768 . 


21 


450 . 


3000. 


-10. 


16.83 


5484. 


16.87 


5481. 


22 


450 . 


3000 . 


-10. 


17.96 


3299. 


17 . 98 


3296 . 


23 


4 50. 


3000 . 


-10. 


10.67 


10940. 


10 . 69 


10935 . 
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FORTRAN AND PL/M RESULTS (continued) 



IDNO 


VKTS 


ALT 


DEG 


FORTRAN 


PL/M 


TOF 


RANGE 


TOF 


RANGE 


24 


450. 


3000. 


-10. 


7.82 


12814. 


7.85 


12815 . 


25 


450. 


3000. 


- 10. 


31.70 


1657. 


31.70 


1656. 


26 


450. 


3000. 


-10. 


36.29 


1384 . 


36.30 


1382. 


27 


4 50 . 


3000 . 


-10 . 


41.58 


1213. 


41.58 


1210. 


28 


450. 


3000. 


-10. 


20.56 


2444. 


20.57 


2442 . 
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EXIT 2: T EMP = ENT (LOOK I ) ; TEMP1=EN T (L00KI+1) ; TEMP2=ENT (LCOKI+2) 

/* CHANGE THE SIGN OF TEMP */ 

T EMP2=TEMP2 XCR 80H; 

/* ENTRY - ENT (LOOK) */ 

CALL ADD( .TEMP, .ENTRY,. TEMP); 
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